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We consider Brans-Dicke type nonminimally coupled scalar field as a candidate for dark energy. In 
the conformally transformed Einstein's frame, our model is similar to coupled quintessence model. 
In such models, we consider potentials for the scalar field which satisfy the slow-roll conditions: 
[(l/V)(dV/d(j))] 2 « 1 and (l/V)(d 2 V/d<fi 2 ) « 1. For such potentials, we show that the equation 
of state for the scalar field can be described by a universal behaviour, provided the scalar field 
rolls only in the flat part of the potentials where the slow-roll conditions are satisfied. Our work 
generalizes the previous work by Scherrer and Sen [3] for minimally coupled scalar field case. We 
have also studied the observational constraints on the model parameters considering the Supernova 
and BAO observational data. 



I. INTRODUCTION : 

Over the last few years, there are growing evidences 
in favour of the scenario that the universe at present is 
expanding with an acceleration. The supernovae type 
la observations the cosmic microwave background 
radiation (CMBR) probes [2| suggest this acceleration 
very strongly. This is confirmed by the very recent 
WMAP data as well. These observations indicate 
that around 70% of the total energy density of the uni- 
verse is in the form of an exotic, negative pressure com- 
ponent dubbed as dark energy (See reference 0] for a 
recent review) . The ratio of pressure to density for the 
dark energy component, referred to as the equation of 
state parameter aide, is given by 
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Recent observations suggest that io d e is close to — 1. 
And assuming the equation of state parameter to be 
a constant, the approximate bound on ujd e is —1.1 < 
Wde < —0.9 (See references d, @]). With future obser- 
vational data, if ojde continues to be around —1, then the 
most reasonable choice for dark energy will be a cosmo- 
logical constant. However, if the observations predict a 
value close to —1, but not exactly equal to —1, dynam- 
ical dark energy models, like those which are driven by 
a scalar field can be an interesting option to look for. 
Such models usually depict the evolution of uide with 
time. This fact has recently been explored by Scherrer 
and Sen 0, H| , who examined the quintessence models 
as well as phantom dark energy models where the scalar 
field potentials satisfy the slow-roll conditions : 
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The advantage of these models is that they provide a 
natural mechanism to produce an equation of state pa- 
rameter LUde slightly less than —1 at present. However 
the disadvantage of most of these scalar field models is 
that the quintessence potentials do not have any sound 
field theoretic background explaning their genesis. So, 
it might appear more appealing to employ a scalar field 
which is already there in the realm of the theory. This 
is where non-minimally coupled scalar field models be- 
come crucial where the scalar field is not put in by hand, 
but is already there in the purview of the theory. Al- 
though a long range scalar field as a gravitational field 
was first introduced by Jordan Q , this kind of non min- 
imally coupled scalar field actually attracted the atten- 
tion of the researcher when Brans and Dicke invoked 
such a field in order to incorporate the Mach's Principle 
in General Theory of Relativity [10]. Brans-Dicke (BD) 
theory, where the scalar field is directly coupled to the 
Ricci scalar, has the merit of producing results which 
can be compared with the corresponding GR results 
against observations [ll[ . A further important virtue of 
BD theory is that it is believed to produce the GR in a 
particular limit (see [l2j for a different result). 

Possibility of having late time acceleration of the 
universe with nonminimally coupled scalar field (a.k.a 
scalar tensor theories) have also been explored in great 
details [l6l |. Scaling attractor solutions with nonmini- 
mally coupled scalar fields have been studied with both 
exponential and power law potentials 17| . Faraoni 
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[181 ] has also studied with a nonminimal coupling term 
(<fr 2 /2)R with different scalar field pot entials for the late 
time acceleration. Bertolo et al. [19], Bcrtolami et al. 
[20l | . Ritis et al. [2l[ have found tracking soultions in 
scalar tensor theories with different type of potentials. 
In another work, Sen and Seshadri [22j have obtained 
suitable scalar field potential to obtain power law accel- 
eration of the universe. Saini et al. [23| and Boisseau et 
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al. [2~ij have reconstructed the potential for a non min- 
imally coupled scalar field from the luminosity-rcdshift 
relation available from the observations. Attempts have 
also been made to obtain an accelerating universe at 
present by introducing a coupling between the normal 
matter field and the Brans-Dicke field in a generalised 
scalar tensor theory [25| . 

In the present work we try to investigate dark en- 
ergy models in non-minimally coupled theories of grav- 
ity, more precisely in BD theory. We write the equations 
in the conformally transformed version of the theory (so 
called Einstein frame) and try to obtain a general form 
of the equation of state parameter ujde- The field equa- 
tions in this version look simpler and are more tractable. 
But now there is a coupling between the matter and the 
scalar field. This is in a way same as the so called Cou- 
pled Quintessence model extensively studied in the liter- 
ature by number of authors 26] . The existence of such 
interaction between the two components can provide 
some ineteresting clue about the mystery of the dark 
energy [27j • The effect of such coupling in different ob- 
servational features like CMB and structure formation 
have also been studied [28]. In the present work,we try 
to obtain an expression for to^e in the BD model in Ein- 
stein frame (which is essentially a coupled quintessence 
model) where the scalar field potential satisfies the slow 
roll conditions @ and ([3]). 

The next section describes the field equations and 
their solutions. Section 3 describes the observational 
constraints on various parameters of the model and the 
last section discusses the results. 
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From equations (O, and §8§ one can easily arrive at 
the energy conservation equation given by 



Pm + "5-pm = , 



(9) 



which is not an independent equation but follows from 
Bianchi identity. 

Let us now effect a conformal transformation 



9nv = e^g^ 



(10) 



where C = 2 " 9 +3 and Imp — The conformal co- 

moving time and the scale factor are now related to the 
original ones through 

±_ - £_ _ 

dt = e 2 Si dt and a = e 2 ^ a . 
The relevant field equations in the new frame look like 

1 ., 
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II. FIELD EQUATIONS AND SOLUTIONS : 

The effective action for Brans-Dicke (BD) theory 
along with a self interacting potential is given by 



S = 



to 
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+ L, 



(4) 

where R is the Ricci scalar, ip is the BD scalar field 
which is nonminimally coupled to the gravity sector, 
a; is a dimensionless parameter called the Brans-Dicke 
constant, V(ip) is the potential for the Brans-Dicke 
scalar field and L m represents the matter Lagrangian. 
Here we have chosen the unit 8ttG = c = 1 . 
For a spatially flat FRW universe, the line element is 
given by 



ds 2 =dt 2 -a 2 (t) [dr 2 +r 2 d6 2 
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Variation of action Q with respect to the metric com- 
ponents and the scalar field respectively yield the Ein- 
stein field equations and the evolution equation for the 
scalar field as 
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where p m = e 7? p m and V = e szV. An overbar 
indicates quantities in the new frame and from now on- 
wards an overdot will indicate differentiation with re- 
spect to the transformed time coordinate i. The field 
equations in the new frame are similar to that for cou- 
pled quintessence model previously studied [2q . 
In the new frame, the energy conservation equation also 
gets modified as 
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where W = y | . Equation ([TJ| shows that the non- 
minimal coupling of the scalar field in the Brans-Dicke 
theory leads to an explicit energy transfer between the 
scalar field and the fluid in the Einstein frame. Fol- 
lowing references @, H§], equations (fTTj) - (fT3j) can be 
written in the form of a plane-autonomous system by 
introducing the variables x, y and A defined by 
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A = 
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These give 

£l<f, = x 2 + y 2 , 
and the equation of state parameter as 
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Following the same method as ref. 0, equations (fTTj) 
and (|13p in a universe containing matter and scalar field 
(as we are mainly interested in the late time behaviour 
of the universe, we neglect the radiation component), 
become 
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x 1 = -3ac + A J -y 2 + -x(l + x 2 - y 2 ) + W(l - x 2 



V = -AW ^xy + '^y(l + x 2 - y 2 ) 



A' = -V6A 2 (r - l)x , 
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where V = V^r/ {^^j and a prime indicates differ- 
entiation with respect to In a. Since we are interested 
in models where lu^ is very close to — 1, it is convenient 
to express in terms of 7, which being very close to zero, 
one can expand various quantities to lowest order in 7. 
Now using equations IT8|) and (fl9|) . we change the de- 
pendent variables x and y to the observable quantities 
fi^ and 7 and obtain 
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It is interesting to note that equations (f2"3"|) , ([2^)1 and 
(|25p are similar to those in ref Q, except that now 
we have an additional interaction term involving W 
arising due to the non-minimal coupling. If we switch 
off this interaction, i.e, W — > 0, this indicates that 
the Brans-Dicke parameter u> — > 00 and we get back 
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FIG. 1: Plot of A as a function of scale factor a. W has 
been choosen to be 0.05. The lower set is for Ao = 0.1 
and the upper set is for Ao = 0.3. For each set the upper 
line is V(<j>) = <fi 2 and the lower line is for V(cf>) — 4>~ 2 ■ 



the GR limit. This is consistent with the notion that 
Brans-Dicke theory reduces to GR in the infinite u> 
limit. 



From equations (|2"3")l and 



one arrives at 



drf_ -37(2-7) + A(2-7) v /3tTV 

., ' 3(i- 7 )(i-n*)[n* + f + 



3(1 - 7 )(i - fy) [n 4 + f 7271^(1 + 7 )] ■ 

(26) 

We should also emphasis that this equation is valid only 
when dQ^/dt ^ (See [3C| for models where oscil- 
lates in time). 

At this point we make two assumptions - the first one 
being 7 << 1, which corresponds to uj^ very close to —1 
as discussed earlier. The second assumption we make is 
that the scalar field starts with an initial value in a po- 
tential which is nearly flat and the slow-roll conditions 
©, © are satisfied. 

Assuming that these conditions are satisfied, one can 
use (|25|) to show that A is approximately constant dur- 
ing the evolution of the scalar field for different poten- 
tials, i.e A = Ao where Ao is some initial value corre- 
sponding to the inital scalar field value (f>o- In figure 1, 



FIG. 2: Plot of W0 vs. f2 for W = 0.05 and A = 0.1. 
Different lines are for V{<j>) = <fi 2 , V(4>) — (f)~ 2 and 
our analytical expression (f2"5| . Thev are all practically 
indistinguishable . 



we have shown that this is indeed true. For a particular 
value of W, W — 0.05, we have shown that the varia- 
tion of A during the entire evolution is very small from 
its initial value Ao for two completely different poten- 
tials e.g V{4>) = (j) 2 and V(<f>) = <f>~ 2 : for Ao = 0.1, A is 
practically constant, while for Ao = 0.3, the variation is 
^ < 2.5% from a = 10~ 3 till today. 

So, we replace A by Ao in equation 126[) and retain 
terms only to the lowest order in 7. Also we retain terms 
upto lowest order of W. Remember that W is inversely 
related with the BD parameter oj which has to be large, 
resulting smaller value for W. This also ensures that 
the departure from the GR result is minimal. Using 
this, we obtain following equation: 
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With W — 0, one can get back to the solution 7(f^) 
as obtained earlier by Scherrer and Sen 0, [!]• With 
W one can also integrate the above equation to get the 
corresponding 7(f2^). In order to do that one need to 
have some boundary condition to determine the inte- 
gration constant arising due to integration of (|27| . The 
ideal condition should be 7 = for = 0. That is to 
say we start at very early time, when the dark energy 
contribution is zero and put the scalar field at the flat 
part of the potential. We can not take the limit = 
initially as term involving W in the expression for 7 af- 
ter integration, blows up. Instead we assume that for 
some very small initial value of f^, say Qi, 7, = 0. We 
should stress that f)j can be very very small but not 
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FIG. 3: Plot of w vs. for W = 0.05. The three 
sets are for Ao = 0.1,0.2,0.3 from bottom to top. For 
each set, top line is for V(</>) = 2 , the middle line is 
our analytical expression (|2"5)) . and the bottom line is 
for V(<j>) = (j)- 2 . 



exactly zero. With this initial condition, one can get 
the equation of state: 
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This is the main 
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result we obtain in this study. This gives an analyti- 
cal expression for the equation of state parameter for 
the dark energy as a function of its density parameter 
for a nonominimally coupled scalar field model (which 
is essentially a coupled quintessence model in Einstein 
frame) . This generalizes the previous result obtained by 
Scherrer and Sen [3, Q for a corresponding minimally 
coupled scalar field model. The expression has an extra 
parameter W which is related to the nonminimal cou- 
pling (or in other words related with the energy trans- 
fer between the matter and the dark energy in Einstein 
frame). The parameters Ao and Hi can be related with 
the present day value of the equation of state parame- 
ter w^o of the scalar field and the corresponding value of 
its density parameter fi^o- Hence effectively the model 
has three parameters W, lo^q and Q^q. Next we want to 
show that our analytical expression (|28[) for the equa- 
tion of state parameter uj^ indeed represents the true 
behaviour for lo^ for variety of potentials V{4>) subject 
to the condition that all of them have a very flat part 
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FIG. 4: Plot of uj^ vs. £1$ for A = 0.3. The three sets 
are for W = 0.05, 0.1, 0.2 from bottom to top. For each 
set, top line is for V(tf>) = <fi 2 , the middle line is our 
analytical expression (|28p. and the bottom line is for 
V(cj>) = 0- 2 . 

over which the scalar field rolls, i.e, Ao is small. One 
peculiar feature of the nonminimal coupling (or coupled 
quintessence) is that even though we start with a scalar 
field which is initially confined to the flat part of the 
potential leading to oj^i — — 1, it starts increasing to a 
higher positive value initially before coming back to a 
negative value near — 1 . And this is not only for our ap- 
proximate analytical expressions for given in (|28p , it 
is also true when we solve numerically the exact equa- 
tions (E2), and (UHl). And this is solely due to the 
interaction between the matter and dark energy. This 
interaction term is related with matter energy density 
(see equation (fT3"]0 and as initially it is large, the in- 
teraction is also very large. This actually violates our 
assumption of 7 << 1 initially. But we should emphasis 
that this does not create any problem for the viable cos- 
mology as this peculiar behaviour happens only when 
has negligible contribution to the total energy den- 
sity of the universe, and does not contribute in any way 
to the cosmological evolution. Once O^, starts con- 
tributing substantially to the total energy density, the 
equation of state wj, reaches value close to cosmologi- 
cal constant which is what we expect. In Figure 2, we 
have shown the behaviour of as a function of f2^ 
for negligibly small values of fi^. We have shown the 
behaviours for analytical expression (|28|) as well as the 
result obtained by numerically solving (|2"3")) - (f2l)|) for po- 
tentials V(4>) = 4> 2 and V(4>) — cf>~ 2 . The plots are 
shown for W = 0.05 and Ao = 0.1 (fii has been choosen 
so as to have fi^, at present to be equal to 0.7). One can 



FIG. 5: Plot of ojtf, vs. scale factor a(t) for Ao = 0.3, 
rj o = 0.7. The three sets are for W = 0.1. The solid 
line is for (|2"9")) . the dashed line for V(<fi) = (j> 2 , and the 
dotdashed line is for V{(j>) = <\>~ 2 ■ As one can see all 
three lines arc indistinguishable. 

see the different behaviours are indistinguishable from 
each other, showing, (i) with small Ao different poten- 
tials have similar evolutions and (ii) our analytical ex- 
pression (|28p represents correctly the exact behaviour 
for small Xq. 

In figure 3 and figure 4, we show the same behaviour 
for larger values of f2^. We have plotted uj^ as 
a function of fta for our analytical result (|2"8|) as well as 
by solving numerically the exact equations (|2"3")l - (f2"5) 
for different potentials, using different values of W and 
Ao. In each case, Hi has been chosen in order to have 
£1^ at present to be equal to 0.7. The figures show that 
for small Ao our analytical expressions nicely represent 
the exact behaviours for completely different potentials. 

Next we have to find the expression for fi^a), in 
order to express lo^ as a function of the scale factor 
a. This is essential for confronting the model with any 
observational result. In this regard, we can use (|23p to 
solve r^(a) to use it subsquently to determine u^(a). 
Assuming both 7 and W to be small, we can neglect 
the terms involving those in (|23[) . to get an approximate 
solution for f^: 

Qt=[l + ( 7 }--l)a- 3 }- 1 (29) 

where fi^o is the the present day value of fi^ and we 
normalize the present scale factor to 1. This is same as 
what we get for a ACDM model which is not unexpected 
as we are considering scalar fields which have equation 
of state very close to —1. To show that this form for 
fit/, very closely resembles with what one gets by exactly 
solving the (|2"3"l) - (f2"S")) for small Ao, we show the plots in 
figure 5. Here we plot our analytical expression ([29)) 
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FIG. 6: la and 2a confidence contours in lu^q — fl 
plane. The dashed lines e 
the solid lines are for Snla-I 
text). This is for W = 0.05 



plane. The dashed lines are for Snla data only and plane. The dashed lines are for Snla data only and 
the solid lines are for Snla+Bao data (explained in the the solid lines are for Snla+Bao data (explained in the 



FIG. 7: la and 2a confidence contours in ui^o — 
plane. The dashed lines 
the solid lines are for Snla- 
text). This is for W = 0.1. 



for together with the result one obtains for f2^ by 
solving numerically the exact equations P3"[) - ([2"5"j) . We 
have fixed W = 0.1 and Ao = 0.3 for this purpose and 
assume fi^o = 0.7. One can see that all three lines 
are just indistinguishable, confirming the fact that our 
analytical approximation (|29p for is indeed robust. 
It also shows the contribution of the scalar field to the 
total energy density of the universe is negligibly small 
till a ~ 0.2 and any peculiar behaviour of the equation 
of state ujfj, before this era does not matter in any way. 

So lj28]) together with the (j29f for fi^(a) represents an 
analytical behaviour for the equation of state of a non- 
minimally coupled scalar field with BD type coupling. 
This behaviour is only valid when the scalar field rolls 
over a very flat part of the potential so that equations 
^ and ([3]) are satisfied, irrespective of the form of the 
potential. This is main result of this study and its gen- 
eralizes the similar result obtained earlier by Scherrer 
and Sen 7, 8] for minimally coupled scalar fields. 
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FIG. 8: la and 2a confidence contours in oj^o — H^q 
plane. The dashed lines are for Snla data only and 
the solid lines are for Snla+Bao data (explained in the 
text). This is for W = 0.2. 



III. OBSERVATIONAL CONSTRAINTS 

Next we compare our result for cj^(a) to the Snla 
and Baryon Acoustic Oscillations data. The confidence 
contours are constructed using 60 Essence supernovae, 
57 SNLS (Supernova Legacy Survey) and 45 nearby su- 
pernoave, and the data released by 30 Snla detected by 
HST and classified as Gold sample by Riess et al [31] . 
The combined data set can be found in Ref. [f|. We 
also add the result obtained by Sloan Digital Sky Sur- 
vey regarding the Baryon Acoustic Oscillation (BAO) 
Peak, giving a distance scale at zt, ao = 0.35. In Figure 



6,7 and 8 we have shown the confidence contours for 
W = 0.05, W = 0.1 and W = 0.2. 

Our result shows that cosmological constant is obvi- 
ously very much consistent with the observational re- 
sult. Also, unless one obtains a stringent bound on 
fi m o from structure formation data e.g the shape pa- 
rameter for the matter power spectrum, which will re- 
sult a corresponding bound on f^Oj one can not put a 
strong constraint on uj^o, although BAO data seems to 
provide a strong upper bound for uj^ which is around 
w o < -0.92. 
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IV. CONCLUSION 

We have studied dark energy models with scalar fields 
which are non-minimally coupled to gravity. We as- 
sume the form of the coupling is of BD type. We write 
the system of equations in conformally transformed Ein- 
stein frame where the scalar field becomes coupled with 
the matter sector although minimally coupled with the 
gravity sector. Our model is exactly similar to the cou- 
pled quintessence models considered previously by sev- 
eral authors [26j]. To study such coupled quintessence 
models, we have assumed that the coupling parameter 
W between the matter and dark energy (which is re- 
lated with BD parameter to in Jordan frame) is small. 
This ensures that our model does not deviate much from 
the standard Einstein gravity. 

With a such a system, we consider potentials for the 
scalar fields which satisfy the slow-roll conditions ©, 
(|3|) . In such a scenario, we show that all of them con- 
verge to a universal behaviour given by (J2HJ) and (J2U). 
This work generalizes the previous result by Scherrer 
and Sen obtained for minimally coupled scalar fields, 
to the non-minimally coupled case which can also be 
described as a coupled quintessence model in Einstein 
frame. 

In this model, one has to fine tune the initial con- 



ditions in order to get desired value for the equation 
of state or density parameters at present which is in 
contrary to the so-called tracker model. But here the 
subsequent evolution of the universe is insensitive to the 
form of the potential as long as the slow-roll conditions 
([2]), ([3]) are met. This is in contrast to the tracker mod- 
els where one has to really fine tune the shape of the 
potential in order to achieve the tracker behaviour. 

We have also tested our model with the observational 
data from Snla and BAO peak. Although to get strong 
bound on the equation of state, one needs to have strong 
independent observational constraint on f^Oj still in- 
clusion of BAO data put a strong upper bound on uj^q 
irrespective of the value of W , the coupling parameter. 
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